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Abstract 

Based on the non-local light-ray operator technique, we develop an algorithm 
for the computational calculation of evolution kernels for higher-twist opera¬ 
tors in leading order of the perturbation theory. We compute the evolution 
kernel for the twist-3 operators in the flavour singlet channel. Our result con¬ 
firms the local anomalous dimensions computed by Bukhvostov, Kuraev, and 
Lipatov as well as the non-local evolution kernels of Balitzky and Braun. 

I. INTRODUCTION 

Recently, in deep inelastic scattering (DIS) the first moments of the transverse polar¬ 
ized structure function g 2 {xBj, Q^) have been measured |^. Among the leading twist-2 part 
g 2 ixBj,Q^) contains a non-power suppressed twist-3part, which it is intended to measure 
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at HERMES and the SMC. Although at present the statistic is too low to extract the 


evolution of the twist-2 part of g 2 {,XBj,Q^) with respect to the momentum transfer it 
may be important for future high-precision measurements to know the theoretical prediction 
also for the twist-3 part. This higher-twist contribution is in fact expressed in terms of three 
particle operators sandwiched between the polarized nucleon state; therefore, it possesses no 
simple parton interpretation. Furthermore, also for (higher-twist) operators with hxed spin 
there remains a mixing problem due to the renormalization. Consequently, the evolution 
equation is no longer of the type of a Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation. 

Introducing so-called qausi-partonic operators, Bukhvostov, Kuraev and Lipatov derived 
evolution equations for higher-twist correlation functions and computed the twist-3 evolution 
kernels and the anomalous dimensions for both the non-singlet and the singlet channels 
Their non-singlet result was conhrmed in the axial gauge as well as in the covariant 
gauge by different authors using quite different techniques. The results in [^|[ coincide 
with the original result; however, in [Q the transposed local anomalous dimension matrix 
was obtained. Recently, two more independent calculations []^,^ support the Bukhvostov, 
Kuraev and Lipatov result. 

It should be mentioned that special emphasis was given in (see also [§) to the equa¬ 
tion of motion (EOM) operators. As is expected from the renormalization properties of 
gauge-invariant operators , the result obtained in the off-shell scheme coincides with the 
previous ones in which the mixing with the EOM operators was neglected from the begin¬ 
ning of the calculation. In the minimal subtraction scheme this result is obvious at leading 
order. However, beyond the leading order it is not excluded that a non-proper regularization 
of infrared singularities spoils the general renormalization principles. 

Here we apply an appropriate technique, which is based on the non-local operator product 
expansion introduced by Anikin and Zavialov and the renormalization group equation 
for light-ray operators [|^. For several years Geyer, Robaschik and others applied this 
approach to the calculation of leading twist non-forward evolution equations |l^ relevant 
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for deeply inelastic Compton scattering and, more recently, for the forward twist-3 non¬ 
singlet channel p]. A similar approach was independently developed by Balitzky and Braun 
1^ and among other studies, applied to the one-loop calculation of the non-forward twist-3 
singlet channel. 

Although the singlet result [Q coincides with the local anomalous dimension matrix given 
in Pi (up to obvious misprints in both papers) it is desired to have a third independent 
calculation. On the other hand it is also theoretically interesting to study the evolution of 
power-suppressed contributions. Calculating the hrst few moments of twist-4 operators, a 
hrst step in this direction was done in [Q. It should be stressed that also one-loop order 
calculations are cumbersome and will be more difficult for increasing twist. Thus, it would 
be desired to have an algorithm that allows the use of computer power. Such an algorithm 
will be given in this paper for light-ray operators. 


II. HIGHER-TWIST CONTRIBUTIONS TO TRANSVERSE POLARIZED DIS 


Applying the non-local operator product expansion ||T^ , p!^ ,|5[| the structure functions for 
DIS are factorized in coefficient functions and forward matrix elements of light-ray operators. 
The well-known result for the longitudinal structure function gi reads in leading order 




( 2 , 1 ) 


q=u,d,.. 

The quark distribution functions (containing quark and antiquark) are defined as matrix 
elements of leading twist-2 operators: 

dK 


Aqq{x,Q ) = 




iXK,{xP) 


( 2 . 2 ) 


where the renormalization point fi of the operator is set equal to the momentum transfer Q. 
The twist-2 operators are defined as 


Og\Ki,K2) = ^'y5i^q{K2x), 


(2.3) 
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where x is a light-like vector. Here and in the following we apply for simplicity the light-cone 
gange, i.e. xA = 0. A gange invariant dehnition in a general gange is obtained by inclnding 
a path-ordered link factor. 

The operator prodnct analysis to leading order snggests that the strnctnre fnnction g 2 
can be written in the form of a generalized Wandznra-Wilczek relation : 


92{xBj, Q^) = -92{xBj, Q^)+ f — 92{y, Q^), 

JxBj y 


(2.4) 


SO that the Bnrkhardt-Cottingham snm rnle |^, i.e. Jo dxg 2 {x,Q‘^) = 0, is obvionsly 
fnlhlled. Here g 2 {x,Q^) = gi{x,Q^) + g 2 {x,Q^) is decomposed in the twist-2 part given by 
gi and a remaining, not-power snppressed, twist-3 part g 2 |]^. The twist-3 contribntion 


92{x,Q^) = - 4^9q{x,Q^) 

q=u,d,.. 


(2.6) 


can now be formally expressed as 


Aqq{x,Q^) = - 


dn 


sp(p,s - {k 


-K 


P,S) e 


iXK{xP) 


X J 2ti{xP) 

where is the polarization vector and the twist-3 light-ray operator is dehned as 


( 2 . 6 ) 


01p(k2,K2)^ / D'^(t£, KiZ, K2x)l)’q(K2x), (2.7) 


with 


D^{u, Bix, K2 x) = + igA^iJ^iu + B2 u]x), u = 1 — u. 


( 2 . 8 ) 


Unlike Aqq{x,Q‘^) this new twist-3 fnnction Aqq{x,Q‘^) has no simple parton interpre¬ 
tation. In fact, nsing the eqnation of motion {ip — m)'ip = 0 it tnrns ont that the twist-3 
operator (for more details see |Q) can be decomposed in a basis of three-particle operators 
containing also the glnon held and in a mass-dependent two-particle operator. Introdncing 
an appropriate dehnition of a three-particle correlation fnnction the net contribntion to g 2 
in leading order can be written as 


x~g2{x,Q'^) = - Y. 4 

q=u,d,.. 


^■mq{x,Q‘^)+ ^ 


X 


dx Jo 


du uAqq{x, u, 


The three-particle correlation fnnction 


(2.9) 
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Aqq{x,u,Q^)= ^ 


;{P,S\Y^{-ku,ku) + (k ^ -K,)\P,S)\f, 2 ^Q 2 e 


iKx{xP) 


( 2 . 10 ) 


An {xPy 

is even under charge conjugation and depends on the variable m, which gives the relative 
position of the gluon held on the light cone and on the variable x. For 0 < m < 1 the gluon 
held lies between the two quark helds. Because of the support property |x| < Max(l, |2 m— 1|), 
the variable x is then restricted to |x| < 1 and can be interpreted as an ehective momentum 
fraction. The introduced ligth-ray operator ^ 2 ) = 0, ^ 2 ) + 0, Ki) is the 

non-local generalization of the so-called Shuryak-Vainshtein operators ||19||: 


^S^{ki,t,K 2] = ig^qiniX] ^ iF^^{rx)±'g^F'‘^{Tx) x^‘il>g(K2x), 


( 2 . 11 ) 


where is the dual held strength tensor. Furthermore, we introduced a 

two-particle distribution function 

dn Sn 


mq{x,Q^) = 


where the operator 


27r (xP)2 


(P, S\M^{0, k) + {k^ -k)\P, S)\^2^Q2e 


{iK,x{xP)} 


( 2 . 12 ) 


M'’(Ki,K2) = ■tp,{Kix)a“‘‘x2,Y{xD){K2x)il,{K2x), Ua/l = 


(2.13) 


is proportional to the current mass m^. 

Using the dehnition ( |2.10| ) the evolution equation for the twist-3 correlation function can 
be obtained in a straightforward manner from the renormalization group equation (RGE) 
of the non-local ligth-ray operators. From the non-singlet result given in one obtains 
an extended DGLAP equation |TB|: 


d 






5{y - u- uz)-^Pqm{z) 

n \v 


d MQ 2 q(s(Q^) f dz ~ , . b!S f y nii 




z \z 


(2.14) 


dQ2 7 -V / 27J- 

Here, the integration region is determined by both the support of and by the kernel 
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?NS 

2Cf - Ca 


Pqq^iz,u,v) = 


© 1 ( 2 :, u, v)K{z, u, i;)]+ + © 2 ( 2 :, u, v)L{z, u, v) — © 2 ( 2 :, u, v)M{z, u, v) 


+ 


2 


^ 3 ( 2 :, M, v)N{z, u, t>)]+ + P ^ ^ j I + 03 (^, U, v)M{z, M, v) - - z) 


Pqm{z, U, v) = Cf{1 - zf, Pmm{z) = -2CF 1 + 2 : - 


^ 1-z\a)' 

where [A{z, u, i;)]+ = 21 ( 2 :, u,v) — 6{1 — z)S{u — v) f dz' f dv' 21 ( 2 ;', u, v') and 


(2.16) 


© 1 ( 2 ;, M, v) = 6{z)6{u — zv)6{u — zv), © 2 ( 2 :, m, v) = 9{—uvz)9{{l — vz}u)9{{z — u}u), 
© 3 ( 2 ;, u, v) = 9{uz)9{uvz)9{{vz — u}u), 


K{z, u,v) = z 
M{z,u,v) = 


f .. \ 

fu —y u\ 

d{u zv) + 


[v{v-u) 

\V ^ VJ 


vu 


, L( 2 :,M,i;) =-sign(M)—- 2 ;), 


U‘" 


2z{l — zv) ^ sign(n)t> \v . u 




N{z,u,v) = 


-5(1 — z) -\ - 6 {u — zv) 

u[v — u] \ u V 


(2.16) 


Note that, due to the evolution, the variable u can no longer be restricted to the region 
0 < M < 1. The moments with respect to y correspond to matrix elements of operators with 
dehnite spin n, so that their evolution equation will be diagonal with respect to n. However, 
for the remaining mixing problem there is no known analytical solution. In two limits for 
a; —> 1 and (only in the non-singlet channel) also in the large limit, the evolution equation 
for g 2 {x, Q^) reduces to an equation of the DGLAP type |2^ . 

In the singlet channel, there appears a mixing problem. The G-even quark operator 


Y^{kj,K2)^ YAk1,K2)^ Y ^5^(^1i0i^2) + 5^(222,O.K,) 


(2.17) 


q=u,d,... q=u,d,... 

will be mixed with the antisymmetric gluon operator: 


G"(/2i, K2) = gr^<^x^x^x'^F:^{Kix)F^/{0)F^^{K2x), GP{k 2, Ki) = -G^iKi, K2). (2.18) 


III. CALCULATION OF TWIST-3 EVOLUTION KERNELS 

The relevant Feynman diagrams for the calculation of the evolution kernel to leading or¬ 
der are shown in Fig. 1. In addition to the one-particle irreducible diagrams, the application 
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of the equation of motion can be represented by reducible diagrams [Q. In the following 
we describe shortly the computational calculation of the twist-3 evolution kernels in the 
light-cone gauge. The applied method can be easily modihed for a covariant gauge and can 
be extended to twist higher than 3. 

Because of the light-cone gauge all operator vertices 

a(Ki,« 2 ) = with v = {^S,-S,G} (3.1) 


have only a pure exponential ^-dependence. Here includes the Lorentz, flavour, and 
colour structure of the operator. However, the gluon propagator contains a spurious pole: 


.g^'^ — {k^x‘' -|- k^x^) /xk 
k"^ + ie ’ 


(3.2) 


which will be regularized by the Leibbrandt-Mandelstam prescription 


kx^ 


(3.3) 


kx {kx){kx*)+ie 
with the dual light-like vector x*, i.e. x‘^ = x*"^ = 0 and xx* ^ 0. 

The Dirac algebra can be performed with one of the usual high-energy programs (such 
as FORM, FeynCalc or Tracer). The integration over the loop momentum and the hnal 
simplihcation of the result were done with a rule-based program written in Mathematica. 
After cancellation of common scalar products between numerator and denominator, as well 
as the decomposition oil/xk...l/x{k+pi), each diagram is expressed in terms of the following 
dimensional regularized tensor integrals: 

/«-«(a bcaK)= f (3 4) 

With the help of the -I— prescription 

the integral (a, b, c, a; k) can be decomposed in 

a—1 

(q, Q,. Q,. ^ b, c, a — P; K = 0), (3.6) 

/ 3=0 
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where 1+^ {a,b, c, a; k) contains no more spnrious poles. The remaining integrals with 

K = 0 give only a trivial contribntion to the evolntion kernel. Snch integrals are common 
in light-cone gauge and can be found in |2^. Since the integrand of (a, 6, c, a; k) is 

analytic in xk, and vanishes, the momentum integration results in 

/ 1 \ i(—1)“"*"^"*"^ r(a-|-6-|-c — j — 2-f e) fy . ^a-i 


b-1 c-il 9 

xy z -- 


1 d 


2db 


Ml 


2db 


— {b"^ - D{y,z)} 
'n 


— {a+b+c—j—2+e) 


—oiiKixb) 




b^^=B^ (y^z) 5 


where B‘‘(y, z) = j/pf + rpj, D(y, z) = ypl + zpl. 


(3.7) 


Performing the derivation with respect to 6^. and taking into account only the ultra-violet 
(UV) divergent part provide the result for the corresponding Feynman diagram. Unfortu¬ 
nately, the output is cumbersome and must be simplihed. In a hrst step all exponentials 
appearing in the result of a given Feynman diagram will be transformed in a unique form 
in such a way that the measure of the parameter integrals remains unchanged, for instance: 


^-i{[K.\{l-y)+K2y\xpi + [K2{z+y)+Ki{l-y-z)\xp2} V, ^ ^-i{[m(i-y)+K2y]xpi + lK2(l-z)+Kiz]xp2} 


. (3.8) 


Finally, such rational functions in the external momenta as 

dy r dz 6(1-y - z) xpi+xp 2 g. 

Jo ^ Jo ^ ^ [yxpi-zxp 2 ]^ ^ ^ 

have to be transformed into momenta-independent ones. Applying partial integration, this 
expression can be written in the desired form 



ry 

r. ^ 1 ., 11 

/ dy 

/ dz 

6{z)—-5{y)- 

h 

Jo 

[ y+ ^+J 


^-i{[K,iy+K,2y\xp-i_ + [K,2Z+K.-i_z\xp2} 


(3.10) 


Summing up the contribution of the diagrams and taking into account the renormaliza¬ 
tion of the quark and gluon helds, the evolution kernels are given by the ^ pole PB| : 

d ^ dx , fy , { f „ Ca 


IX^^Y{ki,k,2) = ^ / dy I dz{ [Cp 


dp? 


y6{z)Y{-Kiy,K2 - Kiy) 


2ti Jo Jo [ \ 2 

2 zY (ki - K2Z, -K2y) [K{y, z)]+Y{niy K2y, K2Z -F kiz) 


+ 
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2z + [N{y, z)]+ - -5{y)5{z) ]Y{ki- K 2 Z, K 2 y) + [N{y, z)]+Y{my, m - mz) 


2 Nfyy 6 {z)Y{Kiy + K2y, Kiy + my) + Cpy'^diz) [M{ki - my) + M{my - ^i)] - 
Nf{m - 1Z2) {^-y - z + Ayz)G{my + ^2?/, mz + mz) + {i - y - z) x 

{G(-(ki - m){y - z), -my - K2y) - G{-mz - mz, {m - «2)(y - z))} |, (3.11) 

d rl py ( 

li^-^G{m, IZ2) = dy dzUl - y - z + 3 yz)[G{my + my, IZ2Z + mz) + 

G{my - k 2, mz - m)] + {l-y - z)[G{m{y - z), m - t^iy) - G{m - mz, -m{y - ^))] + 

g(^ _l_ _/y \ 

[Q{y, 2;)]+-— - 5 {y) 5 {z)) [G{my + K2y, mz + mz )/2 + G{m{y - z), m - kiz)] - 


1 

16 

1 

16 


[2 - 6{y)6{z)][Y{my - m, Kiz - m)-Y{mz - m, f^iy - K 2 )]/m - 
[2 - 12yz5{l -y- z)+ 5{y)5{z)][2Y{my + my, + f^iz)/{m - 1 ^ 2 ) - 


{Y{my - m,Kiz - m) + Y{mz - m,i^iy - «2)}/«i] - (^i ^ m) 

n (i , , a. _ Y , 1 




( 1 -!/)- 


y M^m) 


( 3 . 12 ) 

( 3 . 13 ) 


where we used the following +-prescriptions: 


[A{y,z)]+= A{y,z) - 6{y)6{z) [ dy'[ dz'A{y',z'), ioi A = {K,Q}, 

JO Jo 

[N{y,z)]+= N{y,z) - 6{y)6{z) f dy'f dzN{y',z'), N{y, z) = 5{y - z) —+ 5{z)-, 

Jo Jo y y 

K{y,z) = l + 5{z)^ + 5{y) -, Q{y, z) = ^ (h^)— + ^( 2 /)—) • 
y z 2 \ y z J 

Taking into account the different operator dehnitions, our massless result coincides (up to 
a missing factor Nf in the quark-gluon kernel and an obvious misprint in the gluon-quark 
sector) with the evolution kernels calculated by Balitzky and Braun P, restricted to the 
forward case. The Taylor expansion with respect to m, ^2 provides the local operators and 
their anomalous dimension matrix. If we set the missing upper sum limit in Eq. (55) of p] 
to n, we agree (up to a constant normalization factor in the dehnition of the gluon operator) 
for the massless case with the general analytical expression for the anomalous dimension 
matrix given by Bukhvostov, Kuraev, and Lipatov. Our mass dependent terms coincide 
with their non-singlet result. 
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(d) 



Figure 1: The relevant Feynman diagrams for the calculation of twist-3 singlet evolution kernels 
in the light-cone gauge for: a) the mass operator b) the mixing of the quark operator 
with the gluon operator G‘’, c) the qiiark operator itself and its mixing with the mass operator, 
d) the mixing of the gluon operator with the quark operator, and e) the gluon operator itself. 
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